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' Abstract 

o : 

■ A cosmology is considered driven by a stress-energy tensor consisting of a per- 

fect fluid, an inhomogeneous pressure term (which we call a "tachyonic dust" for 
reasons which will become apparent) and a cosmological constant. The inflation- 
ary, radiation dominated and matter dominated eras are investigated in detail. 



bX). In all three eras, the tachyonic pressure decreases with increasing radius of the 

universe and is thus minimal in the matter dominated era. The gravitational 
effects of the dust, however, may still strongly affect the universe at present 
time. In case the tachyonic pressure is positive, it enhances the overall matter 
5^ \ density and is a candidate for dark matter. In the case where the tachyonic 

pressure is negative, the recent acceleration of the universe can be understood 
without the need for a cosmological constant. The ordinary matter, however, 
has positive energy density at all times. In a later section, the extension to a 
variable cosmological term is investigated and a specific model is put forward 
such that recent acceleration and future re-collapse is possible. 
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1 Introduction 



There are compelling reasons to study a cosmology which is not homogeneous. 
Inhomogeneous models were studied early on by Lemaitre |1, and Tolman j2] and by 
many authors since. Misner for example, postulated a chaotic cosmology in which 
the universe began in a highly irregular state but which becomes regular at late times. 
The models presented here possess exactly this property, which will be realized in a 
later section. That is, in the matter phase, the deviation from FLRW spatial geometry 
is minimal and we show this by calculating the Gaussian curvature of two spheres in 
all phases. The curvature of a two sphere is the same for all values of r in the matter 
domain yielding a three-dimensional space which is isometric to a sphere. Our present 
location may therefore be anywhere in this universe and there is no conflict with 
observational cosmology. The book by Krasihski also contains many inhomogeneous 
models which do not require us to be located at the center of symmetry Some 
more recent studies dealing with inhomogeneous cosmologies include jS], 0, jSji 

0, Uni and nn. 

Inhomogeneous cosmological models are not at odds with astrophysical data. 
It is well known that inhomogeneities in the early universe will generate anisotropics 
in the cosmic microwave background radiation (CMB). Such effects have been studied 
by many groups (|I2], jlHI; [HI; [IS]) using density amplitudes and sizes of inhomo- 
geneities corresponding to those of observed current objects (galactic clusters, the 
Great Attractor and voids). These studies, utilizing a range of reasonable parame- 
ters, have found that temperature fluctuations in the CMB, AT/ (T), ((T) being the 
mean temperature and AT the deviation from the mean) amount to no more than 
about 10~^ — 10~^, which is compatible with observation. Also, arguments to reconcile 
inhomogeneous solutions with cosmological observations may be found in flUJ. The 
inhomogeneity referred to in this paper is a "radial" inhomogeneity compatible with 
spherical symmetry and therefore its effect on the CMB is potentially more difficult 
to detect than the (small) angular deviations. 

In general, at very high energies, our knowledge of the state of the universe is 
highly limited and special assumptions about the matter content and symmetry should 
be relaxed. It therefore seems reasonable to investigate solutions which, at least at 
early times, are less symmetric than the FLRW scenarios. A thorough exposition on 
various inhomogeneous cosmological models may be found in the book by Krasihski 

In section 2 we consider a cosmology consisting of two fluids, a perfect fluid 
(motivated by the successful standard cosmology) and "tachyonic" dust. We use the 
term tachyonic due to the association of this source with space-like vectors in the stress- 
energy tensor. This terminology is also popular in string-theory motivated cosmologies 
commenced by the pioneering works of Mazumdar, Panda, Perez- Lornezana ^H] and 
Sen [HI and studied by many others (see, for example, HE], HI, EDI, EH, El, ES|> 
El, ESI, EHl, EH) EHl EHl and references therein). It should be pointed out that 
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in neither the case presented here nor the string theory motivated case is the source 
acausal as will be pointed out below . 

The tachyonic dust is chosen as a dark matter candidate for several reasons. 
First, it provides one of the simplest extensions to the standard perfect fluid cosmology 
and it is hoped that this model will provide insight into more complex scenarios. 
Second, as will be seen below, the tachyonic dust is a source of pressure or tension 
without energy density and cosmological observations strongly imply that there exists 
a large pressure or tension component in our universe. This pressure also affects the 
overall effective mass of the universe. Multi-fluid models in the context of charged 
black holes in cosmology have been studied in [^0] 

In section El we consider an extension of the model to the case of variable 
cosmological term. We discuss in detail how making this term dynamical affects the 
fate of the universe. 

Finally, this paper utilizes a number of techniques for analyzing global prop- 
erties of the manifold and it is hoped that this will provide a useful reference for the 
mathematical analysis of cosmological models. 

2 Tachyonic dust and perfect fluid universe 

We consider here a model of the universe which contains both a perfect fluid 
and tachyonic dust. This source possesses the desirable properties mentioned in the 
introduction. Namely, the dust contribution is a source of pressure as is required for 
the recent accelerating phase of the universe. A tachyonic dust is the simplest model 
which contributes to pressure and it will be shown that this pressure also makes a 
contribution to the mass of the universe. This field is therefore also a potentially 
interesting candidate for dark matter. 

Aside from spherical symmetry, the sole assumption is that the eigenvalues of 
stress-energy tensor be real. We may therefore write 

T", = [fx{t, r) + p{t, r)] u^u, + r)b% + r)w^w,, (1) 

with 

U^Up = —1, W^Wfj = +1, U^Wp = 0. 

Here fi{t, r), p{t, r) and a{t, r) are the fluid energy density, fluid pressure, and tachyonic 
pressure (or tension) respectively. By comparison of the a(t,r) term in (Q) to the 
stress-energy tensor of regular dust, it can be seen why we choose the term "tachyonic 
dust" to describe this source. Notice that a dust associated with a space-hke vector 
possesses the desirable property in that it yields solely a pressure. It will be shown 
that this tension may produce the observed acceleration of the universe at late times 
PT] , [32] . The source is not acausal as the algebraic structure of is exactly similar 
to that of an anisotropic fluid which is a causal source under minor restrictions and is 
often used in general relativity (see (32], (Hll, [H] and references therein). 
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The time coordinate, t, may be chosen to be coincident with the proper time 
along a fluid streamline (the comoving condition). This gauge, along with spherical 
symmetry, allows a special class of metrics to be written as 



dr'^ 



1 — er^ + eg{r) 
- dt^ + a^{t) da^. 



(2a) 
(2b) 



This form is particularly convenient as one may readily analyze differences between 
models presented here and the standard FLRW models (the e — > limit). Therefore, 
e may be interpreted as the tachyon coupling constant. It is easy to show that ()2all2b|l 
falls in the Tolman-Bondi class of metrics, used extensively in studies of inhomogeneous 
cosmologies. 

Using (PJ and ()2b|l in the Einstein equations with cosmological constant yields: 



^7r/i(t,r) + A =3 



8np{t, r) - A = - 2 
87rQ;(t, r) = 



d(t) ^ 3e e 
d{ty^ 



[r^(r)]' 



d{t) 

W) " [oM J " ^ 2ra2(t)' 



e ^ eg{r} 



(3a) 
(3b) 

2a2(t) [ r2 J ' ^^""^ 
where dots represent partial derivatives with respect to t and primes with respect to 

Enforcing conservation on (0) yields two non-trivial equations: 



er 



9[rj 



r. 



Kt, ^)-* + ^ {3 [fi{t, r) + p{t, r)] + a{t, r)} = 0, 
[p(t,r) + a(t,r)], + ^^M = o. 



(4a) 
(4b) 



Throughout this paper, restrictions a{t) > 0, d(t) ^ and r > are assumed in 
solving the differential equations. In case the tachyon parameter e = 0, one gets back 
the standard FLRW cosmology. 

The orthonormal Riemann components will be useful: 



\-ere 







dit) 
ait) 






1 


a?{t)r 






'a(t)" 
_a(t)_ 



R 



■tete^ 



rd^(t) + er — 



eg'{r) 



2er — eg'{r) 
ra{t)'^ 



(5a) 
(5b) 
(5c) 
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as well as those related by symmetry (hatted indices denote the orthonormal frame). 
The solutions, being local and valid in some domain, need not possess the neighbour- 
hood near r = 0. The singularity at r = will be addressed in a later section. 

In cosmology, two measurable parameters considered are the Hubble parame- 
ter H{t) and the deceleration parameter, q{t). These are: 



Hit) :- 



-Rf 



q{t) := - 



d{t) _ 
a[tj 

a{t)d{t) 



a? it) 



The field equations, (jHaj) . (j3b|l and (jHcj) yield 

A 



\u(+\\^ A e e[rg{r)\ 



6g(t) [H(t)f =87T [fi{t, r) + 3p{t, r)] - 2A 

1 



er 



[H{t)Y[2q{t)-l] =87rp(t,r)-A + 



a2(t) 



2a2(t) 
eg'{r] 



2r 



g[r) 



(6a) 

(6b) 



(7a) 
(7b) 
(7c) 



To study inhomogeneity, the orthonormal Riemann components of the three- 
dimensional sub-space fOaf) are useful: 



R. 



R. 



eg [r) 

2r 
egjr) 



(8a) 
(8b) 



where the tilde is used to denote quantities calculated using the three dimensional 
subspace metric of t = to spatial hyper- surfaces (j2i|). 

Finally, it is useful to define a measure of the inhomogeneity of the spatial 
universe via an inhomogeneity parameter: 



^ ^' R 



r [2er — eg'{r)] 
2 [er^ — eg{r)] 



(9a) 
(9b) 



A homogeneous space is characterized by = 0. Specifically, for the FLRW (e = 0) 
limit, /(r) = 1. 

We next investigate the three major eras of cosmological evolution. 
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2.1 Matter dominated era 



In the universe's recent history, the galaxies which constitute the bulk of the 
ordinary matter have negligible motion relative to the cosmic expansion. Therefore the 
pressure of ordinary matter is approximately zero. Reasonable physics also demands 
that fj,{t, r) > 0. Setting the pressure equal to zero from the equation (j3E|) yields 
(assuming a(t) ^ 0) 



[a{t)d^{t)]' 



Aa\t) 



— Q w 



-C = a constant. 



(10) 



Here, C is the constant of separation. Solving this equation for g{r) one obtains 

e^(r) = (e - C) + e6, (11) 

with h a constant arising from integration. 

The equation for the expansion factor can be analyzed using techniques, many 
of which are well known in cosmology. We include details here for completion. The 
equation, after an integration, may be written as: 



1.2 Mo A 2 C 



(12) 



Here Mq is a constant arising from the integration. In the standard cosmology this 
equation is often compared to total energy conservation and similar equations have 
been studied at least as early as Lemaitre and Eddington [3H]. The terms on 
the left hand side correspond to a kinetic energy, gravitational potential energy and 
vacuum energy respectively. The total energy being constant, (— C/2). The constant 
Mq may therefore be interpreted as an effective mass of the universe and it is of interest 
to investigate how the tachyon affects this constant. 

The equation (fT^ may be used in (jHaj) along with (fTT|) to give the current 
effective mass of the universe: 



Mo = Ai(t, r) 



3 ^ ' 



eb , , 



(13) 



The second term in this equation gives the tachyonic contribution to the effective mass 
of the universe and therefore represents the present mass due to dark matter (which 
is independent of A in this section). 

The fluid and tachyonic energy density and pressures are given by: 



Sufiit, r) 
87ip{t, r) 
87rQ;(t, r) 



A + 3 



eh 



a{t) 



+ 



(14a) 
(14b) 
(14c) 
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Note that the tachyon pressure can be very small today (for large a{t)) although its 
effects through (fT^ can be very large. 

Acceleration in the matter phase may be analyzed by studying the equations 
d7E|l and 

2q-l=^^[C-Aa\t)], (15a) 

H\t) {2q - 1) ~ - A, (15b) 
a{t) 



q =fi{t,r) 



4 a^t) 

— TT- 



6r2a2(t) 3 a^{t)' ^ ' 



3 a?{t)_ 

Note that for positive /x(t, r), q may be negative even with A = 0. This result indicates 
that the tachyonic dust may drive the relatively recent acceleration phase indicated by 
supernova observations [HI], [S2], |HE]- (Recall that positive acceleration corresponds 
to a negative deceleration parameter.) We will discuss in a later section the values 
the parameters in p5cj) must possess for this scenario. However, the emphasis in this 
paper will be on Lambda driven acceleration, the tachyon assuming the role of dark 
matter. 

Acceleration can also be studied by differentiating (fT^ to obtain 

The above "force" equation nicely demonstrates the fact that positive Mq tends to 
produce an attractive force whereas positive A produces a negative or repulsive force. 
The tachyonic effect is inherent in Mq via (fT^ . 

The fate of the universe is governed by the scale factor, a{t). In general, the 
equation for a{t) cannot be solved explicitly. Here we use effective potential techniques 
to study properties of a{t). Figure [T] shows plots of the effective potentials due to the 
matter fields (grey line indicating the function — Mo/a(t)) and the cosmo logical term 
(dashed line indicating the function — g^^T(^) as well as the sum of the two (solid) for 
various signs of Mq and A. 

From the figure it can be seen that for C > situations depicted in figures 
^ (a) and (d) allow solutions which re-collapse even for A = 0. For C < 0, the 
configurations in figures ^ (b) and (d) allow for re-collapse (there are no re-collapse 
solutions for C < if A = 0). In 2 (c) re-collapse is impossible. 

It is of interest to study the geometry of spatial sections generated by this 
solution. As mentioned in the introduction, the space-like hyper- surfaces are not 
surfaces of constant curvature. As well, if we consider the global picture, then the 
parameters discussed can also affect the topology of the universe. Spatial hyper- 
surfaces at t = to possess the line element (equation |2a| except for a scale factor) 

da^ = + dO^ + sin^ 9 d<p\ (17) 

1 — Cr^ + eb 
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a) 



A >0, Mo > 
act) 
3 



A <0, Mo < 





c) 



A >0, Mo < 





Figure 1: Effective potentials in the matter phase. Dashed Hnes denote cosmological po- 
tential, — g^T^Tj) grey lines denote matter potential, — and solid black lines denote net 
effective potential. Re-collapse is possible for scenarios (a), (b) and (d). Parameters to 
produce the graphs are |A| = 0.1 and \Mq\ = 1 although the qualitative picture remains 
unchanged for other values. 



Although ()17p bears a close resemblance to the standard FLRW line element, they are 
not equivalent. The orthonormal Riemann components for ()17|) yield: 



--C- 



(18a) 
(18b) 



and therefore, for e 7^ 0, the three dimensional hyper-surfaces are not of constant 
curvature. For small e deviations are minimal and for e = 0, the hyper-surface is of 
constant curvature C. The inhomogeneity parameter (j9b|) is calculated to be 



Kr) 



_Cf_ 
Cr^ -eh' 



(19) 
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If we wish to treat the solutions as global, then the spatial topology may 
be studied. The two dimensional sub-manifold {9 = vr/2) of the three-metric ()17|) 
possesses line element 

Transforming to the arc-length parameter, I, along a r-coordinate curve, one can 
obtain 



A - Cr^{l) 



(21) 



with A := 1 + eb and [r'( 
solutions are derived: 



r(/) 



A — Cr^{l). Integrating for r{l) > 0, the following 



sm 



2^sinh 
t (/ - lo) 



C{1 - lo) 



for C >0 , A>0 
for C <0 , A>0 
for C = , A > 



(22) 



{lo is a constant arising from integration). It is clear from (PT|) that / is a geodesic 
coordinate. From the periodicity of the sine function, it may be seen that the two con- 
jugate points on the radial geodesic congruences are given by r(/o) = r{lo+'rc /^/C) = 0. 
Thus, one concludes that spatially closed universes correspond only to C > 0. 



2.2 Radiation dominated era 

Here we study the next major phase in the evolution of the universe. The radiation 
dominated phase is characterized by the relativistic fluid equation of state /i(t, r) = 
3p{t,r). Using this along with (jHajl and (|3bj) yields: 



1 [aHt)]- - ha^t) = + ^ lr"'!,{r)]' = -C, 
again C is a separation constant. The equation for g{r) is satisfied by 



(23) 



eg{r) 



,2/5 ■ 



Solving for the scale factor, a{t), one obtains: 

-Ct^ + Kit + K2 



a\t) 



Kie^y^' + K2e""V 3^ + 



Ki sm 



2A 



2j^t 



+ K2 COS 



2j^t 



for A = 
for A > 
+ 1^ for A < . 



(24) 



(25) 



Here, Ki and K2 are arbitrary constants of integration. However, the domain of t and 
the signs of these constants must respect a?{t) > 0. 
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The densities and pressures are given by 

2 



A + 3 



a(t) 



+ 



a2(t) 



C - 



eh 



5^12/5 



57r 



87rp(t,r) =— /i(t,r) 



87ra(t, r) 



6e6 



5^12/5^2 (i)- 

The Hubble parameter is calculated to be: 



(26a) 
(26b) 
(26c) 



H{t) 



1 Ki-2Ct 

2 -Ct2+Klt+K2 



Ki cxp(^2 



I-3C/A 



Kl COS 2 



-K2 sin( 2-y/^tj 



Kl sm(2A/^t)+K2 cos(2A/^t)+3C/A 



for A = 
for A > 

for A < . 



(27) 



The deceleration parameter is provided by: 

(Ki-2Cif 



6reiK2A+3c[Ki exp(2-y|t)+re2Cxp(-2A/|t)j+AKf cxp(4|t)+AK|exp(-4|t) 









2 


A 


ftie V ■i —K2e * 















K2+K2+|^Kisin(2^/^t)+K2Cos(2Y^^t)j +(^) [ki sm(2y^t)+K2Cos(2^/^t)] 



[ki sin (2-^/^4) +K2 003(21/^4) 



for A = 
for A > 

for A < 



Finally, the inhomogeneity parameter of the spatial hyper-surfaces is given by: 



/(r) 



5 Cr 5 + e6 



5 f Cr 5 — eh 



(28) 



The presence of the tachyon affects the spatial geometry. Here spatial geom- 
etry is again studied via the arc-length parameter /. The geodesic equation along an 
r-coordinate curve yields 



dr{l) 
dl 

V{r) ■ 



+ 2V{r{l)) = 1 



2 2r2/5 ' 



(29a) 
(29b) 



One may analyze ()29aj) via similar "effective potential" techniques as in the dynam- 
ics. For positive C, r{l) is bounded regardless of the sign of the tachyonic potential, 
eh (spatially closed universe). For negative C, all allowed solutions are unbounded 
(spatially open universe). For C = 0, the spatial universe is also open. 
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2.3 Inflationary era 

We now investigate the inflationary phase. It is generally believed that the uni- 
verse experienced tremendous expansion over a short period of time. There are many 
physical reasons for believing in this scenario and an excellent review may be found 
in jSH]. Some studies of the scalar tachyon's relevance to inflation may be found in 
fn\ . [22], [21]; [21] ■ In the scenario presented here, the tachyon does not play the role 
of the inflaton. However, the inflationary phase provides one possible mechanism for 
the transition from high tachyon concentration to low concentration. 

Inflationary scenarios are generally supported by the equation of state ^{t, r) + 
pit, r) = 0. This linear combination of ([Haf) and (|3E|) yields: 



a'it) [ln|a(t)|]- 
The solution for g{r) is given by 



e — 



[r'9ir)]' = C. 



(30) 



eg{r) = {e-C)r' + - 



As well, the following modes are found for a{t): 

a{t) =(3oe"' 




Here (3q, Pi, P2 and H are constants of integration. 
The Hubble factor is given by 



for C = 0, Po>0, 

for C > 0, pi>0, 

for C < 0, P2>0, 

for C < 0, p2<0, 

for C < 0, P2 = 0. 



Hit) 



H = constant 
v^A tanh [y/p^it - to)] 

^Mcothiyp^it-to)] 
VW\cot\^VW\{t-to) 

[{t-to)-' 



for C = 0, po>0 

for C>0, Pi>0 

for C <0, P2>0 

for C < 0, p2<0 

for C <0, /?2 = 0. 



(31) 

(32a) 
(32b) 

(32c) 

(32d) 
(32e) 



(33) 
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The corresponding deceleration parameter 




for 
for 
for 



C = 
C > 
C < 



^0 >0 

A >0 

/?2<0 

/?2 = 



(34) 



tan2 y^/W\{t - to) 
. 



for 



for 



C < 



C < 



The source terms are: 



87r/i(t, r) 



Snpit, r) 



1 



3C + d'(t) + - -A, 



(35a) 



a2(t) 



87ra(t, r) 



2e6 



(35b) 



r4a2(t) 



from which it can be seen that the tachyon is naturally diluted by the presence of a 
scale factor which increases rapidly. The fluid density and pressures, however, need 
not dilute as their expressions contain terms proportional to a{t)/a{t) which may tend 
to constant (as in ()32a|l and (j32b|) ) or increase (as in ()32d|l ). We demonstrate several 
scenarios next. 

In the figure |21 The graphs on the left represent the scenario with C > 
("closed inflation") whereas the graphs on the right represent the C = scenario 
("flat inflation") at some fixed value of r. Both scenarios are with A = so that 
/i(t, r) represents the energy density of all fields (dominated by the inflaton, with 
minor contributions from other fields) save for the tachyon, whose density is given by 
a(t, r) in graphs a) and b). The space-time coordinates possess units of lO"^^ metres 
here. Note that for an acceptable interval of inflation (approx a few times 10~^^s), we 
have, in the C = scenario, a dramatic decrease in the density of the tachyon field 
but not the necessarily the inflaton field. In this scenario, inflation must terminate 
by the standard phase transition of the inflaton field. At the end of inflation, the 
tachyon density is much smaller than the densities of the other matter which will be 
dominated by radiation leading to the radiation era. In the C > scenario, both 
a(t,r) and n{t,r) vary with time although n{t,r) (initially primarily the inflaton) 
approaches a constant value while a(t,r) decreases as a^'^it) (this is not obvious from 
graph c, however it can easily be seen, by examining the analytic expressions for a{t, r) 
and iJ^{t,r) with C > 0, that f^(t,r) possesses a term which does not decay with time 
whereas a(t,r) does not possess such a term). It is a simple matter to show that 
parameters exist to produce an increase in the expansion factor by many orders of 
magnitude. The figures |21 show this although their time axes have been truncated to 
show the behaviour of a{t) more clearly. 

As inflation progresses, both models yield a tachyon density whose value de- 
creases to a smaller value than n{t,r). This value can be made small enough as not 
to intertfere with the physical processes that must have occured during the radiation 
dominated era. 
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a 

tachyonic density 



b 

tachyonic density 




Figure 2: Inflationary scenarios: graphs on the left represent a C > model whereas graphs 
on the right represent a C = model. Space-time coordinates are measured in 10^^'^m and 
densities are scaled accordingly. Graphs a) and b) represent the evolution of a{t,r), graphs 
c) and d) the evolution of fi{t, r) and graphs e) and f) the increase in the scale factor (see 
text for explanation). 



The spatial geometry is again studied using the arc-length parameter, /, as in 
the matter dominated era. In this case 



r{l) 



^|(l + 4Ce6)^/^sin 2VC {I - k) + l} for C > 

-^{(4|C|e6-l)^/^sinh 2^^(/-/o)] -l}^^^ for C<0 (36) 
^ sj{l-lof-eh 



for C = . 



Here we see that, from periodicity of the sine function, C > again can yield the 
closed spatial universe. 

Finally, the inhomogeneity parameter in this phase is 



/(r) = 



Cr^ + eh 
Cr^-eh' 



(37) 
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3 An extension to variable Lambda cosmology 



Recent experiments suggest that the universe is presently in an accelerating phase. 
If one accepts that the net mass of the universe is positive, then the present accelera- 
tion can be explained by the figure 2a alone. Thus, the choice A > must be made. 
In case A > is a constant, re-collapse is incompatible with acceleration. There- 
fore, we consider the generalization of the previous sections to the variable A{t) case. 
This scenario has relevance in light of recent models (mainly based on supergravity 
considerations) which predict that the dark energy decreases and that the universe 
re-collapses within a time comparable to the present age of the universe (see [1^ and 
references therein). 

Time dependent fields with equation of state p{t) ~ —fi{t) have been employed 
in the literature to explain certain evolutionary periods requiring positive acceleration. 
There are also compelling reasons from particle physics for treating the cosmological 
term as a dynamic quantity (see HT], [321; references therein). 

The field equations ()3ap. ()3b|). (jHcj) formally remain the same with the excep- 
tion that A = A(t). However, the conservation equation (jlaj) needs to be augmented 
by an additional term. The definitions of the matter, radiation and inflationary phases 
are retained exactly as before. Therefore, the equations for g{r) in all three phases 
remain intact. 

The solutions for g{r) can be summarized as: 




(38) 



with 




for the matter phase 
for the radiation phase 
for the inflationary phase . 



(39) 



The three-geometries are specified as 



da' 



+ d9' + sin' 9 d(j)', 



(40) 



ueb 



(41a) 



eh 



(41b) 



from which one obtains 




(42) 
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The field equations are 



87r/i(t, r) + A(t) 
87rp(t,r) - A(t) 
8na{t, r) 
with conservation laws: 

k{t) 



3 [a\t) + C] 



iu-Deh 



„u+2 



{u + 2) eb 
2r^+2a2(t)' 



2a(t)a(t) + a^(t) + C 



ueb 



(43a) 
(43b) 
(43c) 



/i(t, r),t + ^ + ^^(t) {3 [/i(t, r) + p(t, r)] + a{t, r)} = 0, 



[p(t,r) + a(t,r)]^ + 



(z/ + 2)e6 
^7rr^+%2(t) 



0. 



(44a) 
(44b) 



The dynamical quantities are given by 



d{t) 



[l2{t, r) + a{t, r) + 3p{t, r)] a{t) } + -A{t)a{t) 



1 



(45) 



The Hubble parameter and the deceleration parameter are furnished as 

A(t) _ _C_ _ {u-l)eb 



H^{t) [2q{t) - 1] =8np{t,r) - A(t) 



a2(t) 



C + 



ueb 



(46a) 
(46b) 



As far as experimental evidences are concerned, the matter domain is the most 
relevant. Therefore, the maximum information possible will be elicited from the field 
equations for that domain. Setting p(t, r) = u = and integrating ()43b|) yields the 
"energy" conservation equation: 



-a^(t) — 

2 ^ ' 6a t 



A(r) d[a%T)] 



t2 



Mo 
a{t) 



C 
" 2"' 



(47) 



Again Mq arises from the integration and represents the total effective mass of the 
universe. Furthermore, t2 is another constant representing the beginning of the matter 
phase. In terms of the matter fields, the mass is: 



Mo =/i(t, r) 



-na^{t) 



a3(t) 
6 



A{T)d[ 



to 



-7ra^(t2 



(48a) 
(48b) 
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Here, a possible interpretation is that the first term represents the total mass of 
observed matter ("normal" matter), the second term the tachyonic contribution to 
the dark matter (non-baryonic mass for pressure, "dark energy" for tension) and the 
third gives rise to potential "dark energy" responsible for acceleration. 
The acceleration is provided by (j45|) and (j47|) as: 

= - Jfe + '^HtHt) - ^ Jl A(r) d [a3(r)] 

= -{l^^^^^^^i^^ (49) 

In case ii{t,r) > 0, efe > and A(t) > 0, the above terms on the right hand side 
produce a combination of both attractive and repulsive forces. 

The Hubble parameter in the matter domain is provided by: 

H\t) J-^ r) + a(t, r)] + ^ - (50a) 



2Mn C 



a3(t) a? it) 3a3(t) 
The deceleration parameter in this domain is 

H\t)[2q{t)-l] = ^^-m. (51) 

It is clear from (j51|) that q can be positive, negative or zero depending on 
the values of C and A(t). A specific model will be proposed which accommodates 
a spatially closed, re-collapsing universe with an accelerating period in the matter 
domain. In the cosmology presented here, this may be realized by setting A(t2) > 0, 
C > and e6 > 0. Observations indicate that C has a value very close to zero. 
A C < universe has cubic divergent volume at all times save the origin when the 
volume is zero. If, however, C is extremely small yet positive, one has finite large 
volume in the matter domain without contradicting observations. 

The time periods for inflation, radiation and matter dominated eras are [e, ti], 
[ti,t2] and \t2,T /2] respectively. The time T/2 indicates the initiation of re-collapse 
and thus represents the half-life of the universe. Of course, the boundaries separat- 
ing the domains are not sharp as we have indicated and therefore the above simply 
represents a rough guideline. 

A possible evolutionary scenario is depicted in figure El Here we plot both the 
scale factor a{t) and cosmological term A(t) as a function of cosmic time. The scale 
factor increases greatly during the inflationary phase (in the model presented in the 
figure, the inflation is driven by some matter field, not the cosmological term). This is 
followed by a decelerating phase and, near the present time, a period of acceleration 
follows. This scenario is based on the tachyonic positive pressure model and therefore 
this acceleration is A driven. To allow for re-collapse, the cosmological term decays 
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(starting at t = t^) so that d{t) becomes negative causing deceleration and eventual 
re-collapse. The figure is symmetric about T/2. Furthermore, one may have a cyclic 
universe where the scenario repeats after the "big crunch" . 




Figure 3: A possible scenario for the evolution of the universe. The present time is denoted 
hy t = Iq and the half-life of the universe denoted hy t = T/2. The solid line represents the 
qualitative evolution of the scale factor and the dashed line the cosmological term. 



A suitable A{t) function may be defined by 

m = 



Ao a positive constant 
Ao 
Ao 



for < t < ti 
for ti <t <t2 
for t2<t <h 



(52) 



Ao-e2(t-t3)^ + e4(t-t4)^ for h<t<T/2 



with e2 > and 64 > 0. 

The expansion factor for this case is given by: 



a{t) 



cosh(v^t) 



3C 



2Ao +'«ie^P 



Aoy- 
3 ^ 



K2 exp 



-2 



for < t < ti 

for ti <t <t2 
for t2 < t < ts , 



(83) 



/(«) - C y 



dx 



, , , . There are enough arbitrary parameters in 

I ^^2^\2Mo-!^ai{t2)\x-^~C 

(j53j) so that a{t) can be joined smoothly in the three phases if one wishes to enforce 
sharp boundaries between the phases. 

The function a{t) satisfies the formidable integro-differential equation 



r..,xi2 Ao 2.,^ 2Mo - K^a^t) 







(54) 
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in the interval < t < T/2. The above equation is too difficult to solve analytically 
at this stage. 

The spatial geometry for C > is governed by r(/) as 

da^ = df + (/) [de^ + sin^ 6 dcj)^) , (55) 

1/2 




I (1 + ACehf'^ sin \^VC {I - lo)] + l} 



for inflation 



F ^ (/ — Iq) for radiation (56) 



y^sin [v^(/-/o) 



for matter . 



dr 



, — ■ By previous discussions, in all phases the physical uni- 

verses are closed. Moreover, the total volume corresponding to the three dimensional 
spatial sub-manifold in the matter phase is given by 

' [aito)Y I smWC(/-/o) dl = '^^^^ [a^(to)] , (57) 



C 



(Note that in the limit C 0+, the above volume diverges). 

We now wish to address the singularity at r = 0. The two dimensional ge- 
ometries (restricted to = 7r/2) yield: 

dal = — j- + r^d^''. (58) 

These two-dimensional surfaces embedded in a three-dimensional Euclidean space pos- 
sess the following Gaussian curvatures: 

C + ^ >0 for inflation 
K{r, 0) = <( C + > for radiation (59) 

C > for matter . 

In the matter domain, the surface is locally isometric to a spherical surface of radius 
1 / ^/C. However, the original three dimensional spaces in the equation ()41b|) all exhibit 
a singularity at the limit r — >^ 0+. Therefore, some possible global pictures for these 
three dimensional spaces are provided in figure |3] 

In the figurelH one of the angles is suppressed so that latitudinal lines represent 
two-spheres. Two possible scenarios exist; the figures on the left represent the spatial 
manifold for inflation, radiation dominated and matter dominated eras which include 
r = (the left and right points in each figure). The figures on the right have the 
domains in the neighborhood of r = excised. The left and right boundaries are 
therefore identified. Note that as the evolution progresses, the anisotropy of the spatial 
sections diminishes yielding a sphere in the matter domain. This is therefore quite 
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a) b) 




Figure 4: Possible global pictures for spatial sections in the inflationary (a, b), radiation 
(c, d) and matter (e, f) dominated eras. The diagrams on the left include the points corre- 
sponding to r = 0. The diagrams on the right have a neighborhood about r = removed 
and the boundaries identified. Note that in the matter domain the surfaces are isometric to 
spheres yet singularities still exist at r = as indicated by the "hairs" in diagram e. 



compatible with observation. The poles of the sphere, however, are singular or must 
be excised. The singularity appears to be "soft" in that it is of the conical type. Also, 
the equations, being local, are valid in a domain ri < r < r2 which need not include 
r = 0. It is likely that such a singularity would be absent in a quantum theory of 
gravity which would be manifest at high energies. 

4 Compatibility with current observations 

Current observations indicate a universe which is approximately 5% baryonic 
matter, 20% non-baryonic matter and 75% "dark energy" which is responsible for 
the recent acceleration phase. The "directly" measurable quantities in cosmology 
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are fj,(t,r), H and go- Roughly, in the present epoch (and, as we are deahng with an 
inhomogeneous universe, in our neighborhood of the universe) these quantities possess 
the following approximate values: 



H 



2 



/i(to, ro) ^1.6 X lO-^*^ m-2, (60a) 



,d 7.3 X 10^^'^m-^ (60b) 
a{t)\ ' ^ ^ 

go ~ - 0.4 . (60c) 

Here to and ro are the current time and position respectively. 

The deceleration equation p5c|) provides a relationship (using the above pa- 
rameters along with ()14c|) ) between A and a{t,r) (we assume that any time variation 
in A can be ignored): 



a 



(to, ro) = (8.0 X 10^2^ - 6.9 x 10"^^) m'^ . (61) 



If there is no cosmo logical constant, then the second term in this equation indicates 
the approximate value the tachyon tension must possess in our region of the universe 
to drive the observed acceleration. 

If, on the other hand, the tachyon possesses positive pressure (contributing all 
or in part to the non-baryonic dark matter of the universe) then the acceleration is A 
driven. In such a case a{tQ, vq) may take on the following values: 

< a(to, ro) < 4/i(to, ro) ~ 6.4 x lO'^^ m'^. (62) 

The upper limit comes from noting the observational evidence that the dark mat- 
ter contribution is approximately four times the baryonic contribution to the matter 
content. This sets a restriction on the cosmological constant to be of the order 

A = (10-^2) ^-2_ (g3) 

Alternately we may begin the analysis by using equation (jSlj) and solving for 
A (with the parameters quoted above) 



Also, by equation ()50a|) . we may write 

3C 



3H^ = 87r[/i(t,r) + a(t,r)] +Ao 



C 

2.19x10-^2 ^ gvr [/i(t,r) + a(t,r)] + 1.31 X 10-^2 _2_^ (65) 
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(in the last equation has been used.) Isolating the r) +a(t, r) term and using 
in (HHall yields 



3.48 X 10~^^ + 



4fa3(to) 



57r 



C 



1.3 X 10 



-52 



a(to) 



C 



47ra2(to)' 



(66) 



The C/a'^{t) terms represent the present "radius" -squared of the universe. The left 
hand side of (j66|) is an analogue of the present Newtonian density of the universe. The 
above equation is therefore useful in determining the radius of the universe given the 
density or vice-versa. 



5 Concluding remarks 

This paper considers a simple cosmological model consisting of perfect fluid matter 
supplemented with a "tachyonic dust" . The perfect fluid, with positive mass density, 
makes up the ordinary matter as in the standard cosmology. The tachyonic dust 
term is a source of pressure which, interestingly, can increase the effective mass of 
the universe. In this case it could potentially be utilised as a source of dark matter 
although the clustering properties need to be studied. In case the tachyonic dust term 
is a source of tension, it may be responsible for the observed recent acceleration of 
the universe. This model provides the simplest pressure enhancing extension to the 
successful FLRW scenario. At late times, the solution generates a geometry compatible 
with FLRW. 



Acknowledgements 

The authors thank their home institutions for various support during the produc- 
tion of this work. Also, A. DeB. thanks the S.F.U. Mathematics department for kind 
hospitality. A. Das thanks Dr. S. Kloster for useful informal discussions. We thank 
the anonymous referees for helpful suggestions. 



21 



References 

[1] G. Lemaitre, Ann. Soc. Sci. Bruxelles A53 (1933) 51. Engl, trans. Gen. Rel. Grav. 29 
(1997) 641. 

[2] R. C. Tolman, Proc. Nat. Acad. Sci. USA 20 (1934) 169. Reprint Gen. Rel. Grav. 29 
(1997) 935. 

[3] C. Misner, Astrophys. J. 151 (1968) 431. 

[4] A. Krasinski, Inhomogeneous Cosmological Models (Cambridge University Press, Cam- 
bridge, 1997). 

[5] A. Feinstein, J. Ibafiez and R Labraga, J. Math. Phys. 36 (1995) 4962. 

[6] J. Ibaiiez and I. Olasagasti, J. Math. Phys. 37 (1996) 6283. 

[7] J. D. Barrow and K. E. Kunze, gr-qc /9611007 

[8] J. D. Barrow and K. E. Kunze, Phys. Rev. D56 (1997) 741. 

[9] J. Ibafiez and I. Olasagasti, Class. Quant. Grav. 15 (1998) 1937. 

[10] A. Krasinski, Proceedings of the 49th Yamada Conference on black holes and high- 
energy astrophysics Kyoto, Japan, (Universal Academic Press, Tokyo, 1998) 

[11] J. D. Barrow and R. Maartens, Phys. Rev. D59 (1999) 043502. 

[12] J. V. Arnau, M. Fullana, L. Monreal and D. Saez, Astroph. J. 402 (1993) 359. 

[13] D. Saez, J. V. Arnau and M. Fullana, Mon. Not. Roy. Astr. Soc. 263 (1993) 681. 

[14] J. V. Arnau, M. Fullana and D. Saez, Mon. Not. Roy. Astr. Soc. 268 (1994) L17. 

[15] M. Fullana, J. V. Arnau, and D. Saez, Mon. Not. Roy. Astr. Soc. 280 (1996) 1181. 

[16] A. Mazumdar, S. Panda, A. Perez-Lorenzana, Nucl. Phys. B614 (2001) 101. 

[17] A. Sen, J. H. E. P. 0204 (2002) 048. 

[18] G. W. Gibbons, Phys. Let. B537 (2002) 1. 

[19] T. Padmanabhan, Phys. Rev. D66 (2002) 021301. 

[20] A. Feinstein, Phys. Rev. D 66 (2002) 063511. 

[21] A. Frolov, L. Kofman and A. Starobinsky, Phys. Lett. B545 (2002) 8. 

[22] L. Kofman and A. Linde, J. H. E. P. 0207 (2002) 004. 

[23] P. Matlock, R. C. Rashkov, K. S. Viswanathan and Y. Yang, Phys. Rev. D66 (2002) 
026004. 

[24] M. Sami, P. Chingangbam and T. Qureshi, Phys. Rev. D66 (2002) 043530. 

[25] X. Li, D. Liu and J. Hao, \hep-th/020iU^ 

22 



[26] J. S. Bagla, H. K. Jassal and T. Padmanabhan, \astro-ph/0212198\ 

[27] D. Choudhury, D. Ghoshal, D. P. Jatkar and S. Panda, Phys. Let. B544 (2003) 231. 

[28] C. Kim, H. B. Kim and Y. Kim, Phys. Let. B552 (2003) 111. 

[29] G. W. Gibbons, \hep-th/0301iri\ 

[30] A. Das and D. Kay, Can. J. Phys. 66 (1988) 1031. 

[31] A. G. Reiss, et al, Astron. J. 116 (1998) 1009. 

[32] S. Perlmutter, et al, Astrophys. J. 517 (1999) 565. 

[33] M. K. Mak and T. Harko, Proc. Roy. Soc. Lond. A459 (2003) 393. 

[34] B. V. Ivanov, Phys. Rev. D65 (2002) 104011. 

[35] K. Dev and M. Gleiser, Gen. Rel. Grav. 34 (2002) 1793. 

[36] A. G. Reiss, et al, Astrophys. J. 560 (2001) 49. 

[37] G. Lemaitre, Ann. Soc. Sci. Brux. A47 (1927) 49. 

A. S. Eddington, Mon. Not. Roy. Astron. Soc. 90 (1930) 668. 



[39] R. Brandenberger, Proceedings of the International School on Gosmology Kish Island, 
Iran, (Kluwer, Dordrecht, 2000). 

[40] R. Kallosh and A. Linde, J. Gosmol. Astropart. Phys. 02 (2003) 002. 

[41] P. G. Bergmann, Int. J. Theor. Phys. 1 (1968) 25. 

[42] R. V. Wagoner, Phys. Rev. Dl (1970) 3209. 

[43] A. D. Linde, JETP Let. 19 (1974) 183. 

[44] J. M. Overduin and F. I. Cooperstock, Phys. Rev. D58 (1998) 043506. 



23 



